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I. INTRODUCTION 


Heavy-ion experiments at RHIC are now analysing large quantities of data and hope to 
identify and characterise the plasma phase of QCD. As a result it has become imperative to 
put together an accurate and detailed picture of this phase from theoretical analyses of QCD 
at hnite temperature. Present day lattice computations are able to do a large part of this job. 
The reach of theory would certainly be enhanced if perturbation theory is tested accurately 
against these computations and, if found to work well, subsequently used in contexts where 
the lattice is too ponderous to use. 

The gluon sector of the QCD plasma has been extensively analysed. In the quenched 
theory the scaling of has been established to high accuracy The continuum limit of 

the equation of state has also been determined to reasonably high accuracy 11. This shows 
strong departures from the ideal gas. It turns out that hard thermal loop (HTL) resumma¬ 
tion is not sufficient for a description of the lattice results 0, nor is a Borel resummation 
of the perturbative series [^. Other techniques such as an approximately self-consistent 
resummation [^] and dimensional reduction have been applied, and the former has been 
successful in describing the lattice results for T > 2Tc. In addition, the spectrum of screen¬ 
ing masses in the QCD plasma has been determined to good accuracy 0 and, for T > 2Tc, 
found to be in agreement with that determined non-perturbatively from a dimensionally 
reduced model whose couplings are matched perturbatively to the four-dimensional theory 
p. There is also direct evidence for the existence of the gluon as a quasiparticle in the 
plasma for T > l.STc [ITU . 


The quark sector of the theory has not yet been explored in as much detail. Part of 
the reason is that lattice computations with light sea quarks are very expensive. While sea 
quarks are needed in the neighbourhood of the phase transition, where they influence even 
the order of the transition, sufficiently far away their effects are seen to be small 0.0. 
In this paper we report on the physics of valence quarks in the continuum limit of QCD for 
T = 1.5Tc-3Tc, where sea quark effects are marginal. The two aspects we concentrate on 
are quark number susceptibilities and screening correlators and masses. 

Quark number susceptibilities [^, |T^, [T^, |T^, [T^, [T^ now occupy a position of high inter¬ 
est. Due to their connection with fluctuation phenomena and strangeness production 
TSl they are useful inputs to relativistic heavy-ion collisions. They are coefficients of the 
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Taylor expansion of the free energy density of QCD in terms of the chemical potential, and 
hence are useful checks on recent attempts at numerical studies of hnite-density QCD . 
For exactly the same reason they are a test-bed for suggested resummations of the high- 
temperature perturbation series for QCD, which aim to reproduce the free energy density, 
and hence the equation of state, at high temperatures []^. In this paper we report our 
study of these susceptibilities at a large variety of lattice spacings and an extraction of their 
continuum limit. 

Screening correlators and masses are closely related to linear response functions such 
as susceptibilities. They have been explored earlier at finite lattice spacings, 1/4T-1/8T. 
Staggered quarks were used in 4 flavour |2^, 2 flavour ^ and quenched |2^, |2^, |2^ 
simulations. For reasons of symmetry a large degeneracy of screening masses, of the zero 
temperature scalar (S), pseudo-scalar (PS), and certain components of the vector (V) and 


axial-vector (AV) quark bilinears, is expected in the QCD plasma [^. This was not ob¬ 
served; instead lesser degeneracies were— between parity partners S and PS or V and AV. 
Similar results were also obtained with overlap and Wilson quarks in quenched 
QCD, though the splitting between the S/PS and V/AV sectors was seen to be substantially 
smaller. Here we study these correlators on a wide range of lattice spacings and extrapolate 
the results to the continuum. 

In this paper we carefully address the question of the continuum limit of these quantities 
in the quenched QCD plasma, where the valence quarks are used as probes of the dynamics 
of gluons. To this end we perform quenched computations over a range of temperatures and 
a large range of lattice spacings, in order to perform a continuum extrapolation. We test 
the extrapolation by working with two different formulations of quarks on the lattice, with 
different systematics at hnite lattice spacing. One is the normal staggered formulation of 
quarks, the other a Symanzik improved version of staggered quarks, called Naik quarks p8 . 
By this means we identify a unique and verihable continuum limit. 

We would like to distinguish our approach from Symanzik (or loop level) improvement 
of the action m- Such programs aim to approach the continuum limit of the action on 
coarse lattices by adding terms to the action which systematically, order by order in the 
lattice spacing a, remove lattice artifacts up to some sufficiently high order in a. When 
such a scheme is used in conjunction with improved operators, it (ideally) obviates any need 
to go to very fine lattice spacings. In this work we only use an improved operator for the 
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valence quarks; the action is not improved. We reach the continuum limit by simulations 
on a variety of lattice spacings and a smooth extrapolation. 

In the next section we present dehnitions of all the quantities that we measure, and results 
for these quantities for ideal quark gases on a lattice. In Section 3 we present details of the 
simulations and measurement schemes. The following section contains detailed results on the 
lattice with discussion of the extrapolation to the continuum. The continuum extrapolations 
are collected and discussed in the hnal section. 

II. DEFINITIONS 

The partition function of QCD with three quark flavours is 

Z{T,Hu,f^d,l^s)=['DUexp[-S{T)] det M(T, m/,/i/), (1) 

f=u^d,s 

where the temperature T determines the size of the Euclidean time direction, S{T) is the 
gluon part of the action, and the determinant of the Dirac operator, M, contains as pa¬ 
rameters the quark masses, ruf, and the chemical potentials, fxj for flavours /. We also 
dehne the chemical potentials fio = Hu + — hd fis = + f^d — 2/is, 

which correspond to the diagonal flavour SU (3) generators. Note that /io is the usual baryon 
chemical potential and /is is an isovector chemical potential. 

The quark number density for flavour /, Uf, is dehned as the derivative of F = \ogZ/V 
{Z is the partition function and V the spatial volume) with respect to /i/. The quark number 
susceptibilities are the second derivatives 

d^F 

X//,(T,/i„,/id,/is) = (2) 

We will always drop the repeated subscript for the diagonal susceptibilities. There are similar 
susceptibilities, xo, Xs and xs obtained by taking derivatives with respect to /io, /^s and /is 
respectively. Here we determine the susceptibilities at fif = 0. In this limit, all nf{T) = 0 
but the susceptibilities can be non-vanishing. 

The flavour off-diagonal susceptibilities such as 

T 

Xff = y (3) 

are given solely in terms of the expectation values of quark-line disconnected loops. Here 
Mf denotes the Dirac operator for a quark of flavour /, and Mj and M'j the hrst and second 
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derivatives with respect to /i/, taken term by term. The flavour diagonal susceptibilities 
have contributions from both quark-line connected and disconnected pieces— 


Xf 





+ (tr 


(4) 


We draw attention to the fact that the last two terms individually grow rapidly with lattice 
volume, diverging in the infinite volume limit, such that the difference is a finite quantity. 
This cancellation of divergences is the result of a proper treatment of quark chemical po¬ 
tentials on the lattice S' Numerically, the simplest quantity to evaluate is the diagonal 
isovector susceptibility 


T 

y, = — (ti - M~^M'M~^M' 


2V 


(5) 


where the factor half takes care of the isospin carried by each flavour. In addition to this, 
we shall need the baryon number and electric charge susceptibilities 

1 1 
Xo= g (4X3 + Xs + 4Xnd + 4xns) and Xg = g (lOxs + Xs + Xnd - 2Xns) ■ (6) 

Eqs. (^|-^ have been written per flavour of quarks, i.e., the normalisation for the continuum 
ideal gas is Xf/T"^ = 1- 

The screening correlators of local quark-bilinear (“meson”) operators are defined as— 


Cr{z) = Y,{Mj{x,y,z,t)TM\^{x,y,z,t)T) 

x,y,t 


(7) 


where T denotes some Dirac-flavour structure, a and (3 are colour indices and M“^(r) is 
the inverse of the Dirac operator for a point source in the fundamental of colour SU{3) at 
the origin. Since the partition function of eq. (Q) can be expressed as the trace of a spatial 
transfer matrix, the screening correlators decay exponentially, with characteristic length 
scales called screening lengths (inverse of screening masses, fir)- 

The finite temperature symmetries of such a spatial transfer matrix are quite different 
from that of the Hamiltonian and have been worked out in detail [^. In particular, it 
turns out that there are only two independent local correlators. The combination of vector 
operators — Vy and axial vectors AVa, — AV^ are in the representation of the finite 
temperature symmetry group D^. [^. All others, i.e., the scalar S, the pseudoscalar PS, and 
all the components of the V and AV orthogonal to the Bi^, he in the scalar representation. 
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. This clearly implies that the angular momentum J is not a good quantum number for 
screening correlators, and hence states of different J can mix with each other. 


The meson susceptibilities are the zero momentum part of the screening correlator 
and can be written as 


Xr = (trM-iTM-ir) = ^ Z,{T) (pj.) ' , where 


= 




( 8 ) 


|i) denotes the eigenstate of the transfer matrix with eigenvalue A* and the screening masses 
are pp = log(Ai/Ao). In the continuum limit, since M' = 70 , clearly xa = xvo! f^^t, this 


identihcation underlies the computation of |^. The PS susceptibility is also involved in a 
chiral Ward identity— 

{'ip'll)) = mxps- (9) 


At T = 0, chiral symmetry is broken and the sum over states in eq. @ is dominated by the 
lowest term, due to the vanishing Goldstone pion mass. Then this Ward identity shows that 
mlp oc 1/m. However, for T > Tc, chiral symmetry is restored and (ipip) should vanish in 
the chiral limit. At the same time, xps need not be dominated by one state and could 
depend weakly on m, going to a hnite non-zero limit as m —>■ 0. It is interesting to note 
that since the screening masses are equal for all the channels, the differences between 
the susceptibilities can only come from the ZPs in eq. (P). 

For an ideal gas of quarks on the lattice, the Dirac operator at zero chemical potential is 
most simply written in momentum space where it is diagonal— 


Mpg = 5pg [ma + i'y ■ p], where pp = Cio sinp^ G 30 sin 3p^, ( 10 ) 

the momentum spectrum is po = {27r/Nt){n -|- 1/2) with 1 < n < Nt and pi = 2'mr/Ni with 
1 < n < W in all other directions, and pp = sinpp. For staggered quarks, where only a 
nearest neighbour difference appears in the Dirac operator, Gio = 1/2 and G 30 = 0. For 
Naik quarks, on the other hand, Gio = 9/16 and C^q = —1/48. The construction of the 
number densities and conserved currents, and hence the chemical potential on the lattice, is 
complicated for Symanzik improved quarks, but goes through in the usual way 

All the quantities of interest are exactly computable for an ideal quark gas, not only 
in the continuum, but also on the lattice. In terms of 1/D{p) = {ma)‘^ + P ' P-i Pp = 
-(Gio cospo + 3(^30 cos 3po)6op and p/ = {Cio sinpo DGso sin 3po)5o;, we have— 

{'Ip'lp) FFT = £>(?). 

p 


6 














Nt 


FIG. 1: The ideal gas on a lattice with = 3 for staggered quarks (pluses) and Naik quarks 

(crosses), as a function of Nt = 1/aT. 




o Xfft — 




Y. {^(p) \p-p' -f ■ p”] + 2£>^(p) ip ■ p"f] , 


S p 


a/ip'pT = 2 sinh ^ \J(may + sin^(7r/Ait) . 


( 11 ) 


Figure D shows the susceptibilities for the staggered and Naik quarks on a sequence of lattices 
with hxed aspect ratio Ng/Nf = 3. Note that the approach to the continuum limit is quite 
different in the two cases. In free-held theory Xud = 0. The computation of the correlators 
Ct{z) in the ideal gas is also straightforward. Perturbative expansions in the gauge coupling, 
have been performed in the continuum, xs been computed in the hard thermal loop 
(HTL) approximation to order as well as in a skeleton graph resummation [|T^. To the 
best of our knowledge, there has been no such computation of screening correlators, screening 
masses or the remaining susceptibilities. 
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TABLE I: The lattice sizes, Nt x N^, Wilson coupling, (3, and statistics used in this study. The 
statistics is reported as Wep x where Wep is the number of sweeps between measurements 

and Nstat is the number of measurements. 1000 initial sweeps were discarded for thermalisation in 
all cases except the Nt = 14 run, where 7000 sweeps were discarded. For asymmetric spatial sizes 
the long direction is called the z-direction. 

III. METHODS 


As already mentioned, we have investigated quark number susceptibilities using two differ¬ 
ent realisations of lattice quarks. For staggered quarks, the derivative in the Dirac operator 
is discretised using a one-link separated hnite difference operator. This has a discretisation 
error of order a^. The improvement suggested by Naik is to take a specific three-link term 
(see eq. |T^ which cancels the leading term of the discretisation error for staggered Fermions 
up to 0{a^). This possible improvement comes at a fourfold increase in computational cost, 
due to the more complicated structure of the Dirac operator. 

Since the discretisation error in the evaluation of the susceptibilities comes partly from 
the gauge action and partly from the Dirac operator, and we use the Wilson gauge action, 
the part of the O(a^) error coming from this source would remain the same in the two deter¬ 
minations. However, if a large part of the error in the determination of susceptibilities 

using staggered quarks comes from the traces of the Dirac operator in eqs. (M), then the 
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slope in the extrapolation to the continuum limit should be substantially different between 
the two quark formulations. That this is so can be seen already in the computation for the 
ideal quark gas, illustrated in Figure |^. 

In this paper we report results on susceptibilities and screening measured on lattices with 
4 < iVt < 14 for T = l.STc, 2Tc and ST^. The lattices sizes, couplings and statistics used are 
reported in Table Due to statistical errors in the determination of critical couplings various 
other quantities that go into a scale determination, the temperature corresponding to a given 
lattice coupling may be in error by about 5% |^]. The simulations have been performed 


with a Cabbibo-Marinari pseudo-heatbath technique with 3 SU{2) subgroups updated on 
each hit, each sweep consisting actually of 5 sweeps of this algorithm. For susceptibilities we 
have made measurements using both staggered and Naik quarks at m/T^ = 0.1, 0.5 and 1. 
For staggered quarks we have also used m/Tc = 0.03. For screening masses and correlators 
we have used only staggered quarks. 



FIG. 2: The variation of the estimate of on a single configuration with NcGi which values 

as ln(l/ec'G)- The values of ecG are marked at appropriate points on the curve. The size of typical 
errors on the estimates, after averaging over 50 gauge configurations with = 10, is shown as 
the bar. This is insensitive to ecG- 
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The stochastic evaluation of traces depends on the identity 


Nv 




2K 




(12) 


■V i=i 


where X is the identity matrix, Rf is the component a of the vector Ri, which is one of a set 
of N.U uncorrelated vectors with complex components drawn independently from a Gaussian 
ensemble with unit variance, and the star denotes complex conjugation. This leads to the 
stochastic estimators. 


Tr A = Tr AX = ^ y RlARi. 

O AT ^ * * 


2K, 


(Tr A)2 = TtAIAJ 


(13) 


where the vectors used in the stochastic estimators of the two identity matrices, X and J', 
are different. This is done by dividing the set of Ri into disjoint blocks. 

It turns out that the divergence cancellation in eq. is numerically the hardest part 
of the computation. In previous works we had taken a half-lattice version of the staggered 
Dirac operator, which performed extremely well on small lattice volumes but was inefficient 
at divergence cancellation on large volumes , leading to enormous statistical fluctuations. 
We found that the full lattice version of the staggered Dirac operator performs much better 
in this respect, and has no trouble handling our largest lattices. The statistical fluctuations 
are under good control already with W = 10, and one has no need to look for reduced 
variance versions of eq. (0). We found that it is important to work in that sector of the 
quenched theory in which the Wilson-Polyakov line, X, is purely real. This can always be 
achieved by a rotation in the quenched theory. 

The slowest part of the computation is the inversion of the Dirac operator, M, which is 
done by a conjugate gradient iteration. The time taken is linear in the number of iterations, 
Ncg- This in turn depends on the precision required in the solution, which is specified by 
a single number ecc through the requirement that \Mx — rp < NtVecc (here r is one of 
the random vectors and x is the estimate of M~^r). In Figure § we show how the estimate 
of Xs/T'^ changes with Ncg oc log(l/ecG)- If is clear that with ecc = 10“^ the result is 
statistically indistinguishable from that with ecc = 10“®, but the solution is achieved at 
half the CPU cost. We have used ecc = 10“^ in all our computations of x- It turns out that 
the estimation of screening correlators is more sensitive to the stopping criterion, requiring 
ecG = 10-^ 
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The fits and fit errors involved in the continnnm extrapolations are done by a bootstrap 
method. In this, bootstrap samples of the data are extracted and hts are performed on 
statistics obtained from each snch sample. The statistics of the fitted parameters are then 
nsed to obtain non-parametric estimators of the mean and variance. We have varied the 
nnmber of bootstrap samples by two orders of magnitnde, between 10 and 1000, to check 
that the estimators are stable. 

We have nsed staggered qnarks to determine screening correlators and masses. Recall the 
well-known fact that for these, the one link transfer matrix contains complex eigenvalnes and 
leads to oscillations in the screening correlators. However, the two link transfer matrices have 
real eigenvalnes, and the corresponding states can be projected by a simple parity projection. 
Dne to the three-link terms in the Naik Dirac operator, real eigenvalnes are obtained only for 
the fonr-link transfer matrix. The one link transfer matrix contains oscillatory mixtnres of 
varions states from which pure states need to be projected out. We show in the next section 
that we can reach the continuum limit with staggered quarks. Since they lead to simpler 
correlators, we have analysed screening phenomena only with these. We have concentrated 
on degeneracies of the screening correlators as well as the screening masses. The latter 
were extracted from local masses as well as by making two-mass hts to the parity projected 
correlators. Details of these techniques remain the same as in |^. For bias-free extraction 
of screening masses, large separations in the spatial directions are necessary. Consequently, 
we have only used data from lattices with W = 4W for this analysis (see Table |). 

IV. LATTICE RESULTS 

As shown in Figure an ideal gas of staggered quarks has strong lattice artifacts and one 
needs lattices with W — 8-12 in order to get continuum results. Computations in quenched 
QCD, shown in Figure also indicate this. Note specially that for W = 4 and 6, Xs/Xfft 
decreases with W at hxed T, unlike at larger W where it increases. This means that with 
staggered quarks one needs fairly hne lattices for a smooth extrapolation to the continuum. 

Unfortunately, Naik quarks also have their own quirks. An ideal gas of Naik quarks on 
a lattice approaches the continuum smoothly only for W > 6 as shown in Figure |^. The 
interacting theory also exhibits departures from smoothness, as shown in Figure ^ At l.STc 
it might seem that the Nt = 6 and 8 lattices scale, but at 2Tc this scaling is seen to be 
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Nt 

1.5Tc 


2Tc 


3Tc 



Staggered 

Naik 

Staggered 

Naik 

Staggered 

Naik 


0.03 0.1 

0.1 

0.03 0.1 

0.1 

0.03 0.1 

0.1 

4 

1.816 (16) 

1.112 (22) 

1.846 (6) 

1.137 (14) 

1.911 (5) 

1.114 (15) 

6* 

1.312 (59) 

0.957 (36) 

1.428 (5) 1.425 (5) 

1.033 (14) 

1.492 (4) 

0.993 (10) 

8 

1.110 (11) 1.105 (11) 

0.936 (9) 

1.144 (11) 1.142 (11) 

0.952 (8) 

1.191 (19) 1.190 (19) 

0.979 (5) 

10 

1.014 (6) 1.009 (6) 

0.921 (7) 

1.060 (8) 1.058 (8) 

0.948 (4) 

1.066 (7) 1.066 (7) 

0.956 (7) 

12 

14 

0.973 (14) 0.969 (14) 

0.915 (7) 

0.994 (7) 0.992 (7) 

0.985 (10) 

0.941 (9) 

0.949 (2) 

1.027 (6) 1.026 (6) 

0.950 (5) 


TABLE II: Lattice results for with m/Tc = 0.03 and 0.1 for staggered quarks and 0.1 for 

Naik quarks. For Nt = 6 the three temperatures are 1.33Tc, 2Tc and 3.13rc. 



FIG. 3: The temperature dependence of with staggered (left panel) and Naik (right panel) 

quarks of masses 0.03rc and O.lTc respectively, with A) = 4 (circles), 6 (boxes), 8 (pentagons), 10 
(down triangles) and 12 (up triangles). Note the differences in the y-axes in the two cases. The 
bold lines join the centers of the data points. The thin lines show the value of for an ideal 

gas on the same lattices at 3rc. 

accidental. In turn, at 2Tc one sees Nt = 8, 10 and 12 as scaling, but at 3Tc this turns out 
to be accidental. 

Our results for obtained with staggered and Naik quarks for the lightest quark 

masses are collected in Table |1[ In Figure ^ we show the extrapolation of staggered and 
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Naik quark results on xs at T = 2Tc to the continuum limit. For staggered quarks only the 
data for Nt >S hts a simple quadratic extrapolation of the form X 3 (a)/T^ = xs/T^ + 

. In fact, for smaller Nt the extrapolations of Xz/Xfft even gives the wrong sign for s 
. However, for Nt > 8 the extrapolations of and Xs/Xfft give identical results. 

Since the two ratios reach this limit from different directions, their agreement is already a 
good indicator of having the right continuum limit. For Naik quarks at 2Tc, the quadratic 
extrapolation looks reasonable for Nt > 4. This is accidental; at 3Tc or at 1.5Tc a good 
£t is obtained only for Nt > 6. Furthermore, as already mentioned, the extrapolation of 
Xs/Xfft for Naik quarks does not coincide with that of Xs/T^ even at the 3-a level. The 
extrapolation of Xz/T"^ does, however, agree within errors with the continuum extrapolation 
from staggered quarks. 

We reiterate the following point. In this quenched work with the Wilson action we have 
taken staggered and Naik quarks to provide two different operators to measure the same 





FIG. 4: The continuum extrapolation of staggered and Naik quark data for X 3 at 2Tc. The circles 
are the extrapolation of Xs/T^ for staggered quarks, the boxes for XsIxfft for staggered quarks, 
and the pentagons for Xa/T^ for Naik quarks. The bands are 1-a error bands for an extrapolation 
in a? oc 1/Nf. 
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FIG. 5; The parity projected screening correlators evaluated with staggered quarks of mass 0.03Tc 
at T = 2Tc. The degeneracy of PS (circles) and S (boxes) correlators is evident, as is the equality 
of the screening masses for these and the V (pentagons) correlator. The S and PS correlators are 
divided by 3 to bring them into coincidence with the V. The other correlators, which are not 
shown, are also in equally good agreement with these. 



FIG. 6: Local masses for the PS correlator (circles) compared with the result of a two-mass fit. 
The data are for a 12 x 48 x 26^ lattice at 2Tc with m = O.OSTc for staggered quarks. Also shown 
are local masses for the V correlator (boxes) in the representation. 
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Nt 

l.STe 

2T, 

3T, 

4 

4.12 (12) 

4.576 (16) 


6 


5.34 (6) 


8 

5.536 (24) 

5.744 (32) 

5.928 (64) 

10 

5.67 (2) 

6.07 (4) 

6.04 (10) 

12 

5.904 (96) 

6.156 (48) 

6.156 (180) 


TABLE III: local masses for staggered quarks with m = 0.03Tc. 



FIG. 7; Extrapolating the ratio of the measured screening mass, fips/T, and its value in an 
ideal quark gas to the continuum limit. The data are for T = 2Tc and m = O.OSTc for staggered 
quarks. Already at lattice spacings a = 1/lOT and greater the result is compatible with unity. 

physical quantity, namely X 3 . The agreement between the two provides a check that we 
have control over the continuum limit. In this approach we find that equally small lattice 
spacings are needed for the two kinds of quarks, and therefore Naik quarks are fourfold more 
time consuming. The use of improved actions to get continuum results with coarse lattices 
is a different problem. Our observations do not rule out the possibility that improved gauge 
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FIG. 8: The Af+ screening correlators in the PS (circle) and unprojected V (pluses) channels 
compared to their ideal gas counterparts (heavy lines) for staggered quarks. These were determined 
on 12 X 26^ X 48 lattices at T = 2Tc with staggered quarks of m = 0.03Tc. The best fit to the PS 
data is shown (broken line) as is the same fit multiplied by an overall constant to superpose it on 
the data for V, in order to show that these two correlators give the same screening mass. Note 
that the short distance part of the V correlator agrees rather well with the ideal gas expectation, 
unlike the PS correlator which differs by a factor two. 

actions with fat link dynamical Naik qnarks approach the continnnm action more rapidly 
and that with snch actions fat link Naik qnark operators have smoother approach to the 
continnnm already at coarser lattice spacings. 

The off-diagonal snsceptibility, Xud/T"^, is consistent with zero at 95% conhdence level on 
all the lattice sizes, conplings and masses investigated, and at the 67% level on most. The 
extrapolation to the continnnm of Xud/T"^ can be performed by the same fnnctional form 
as for xa/T^, and, nnsnrprisingly, gives a resnlt consistent with zero for both staggered and 
Naik qnarks. For the strange qnark mass region also, similar resnlts are obtained. In view 
of these, we neglect the qnark-line disconnected amplitndes in all the snsceptibilities in eq. 
(|). The continnnm resnlts for Xs/T^, Xo/^^ and Xq/T"^ are discnssed in the next section. 
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We have seen clear degeneracies of the S and PS correlators and the V and AV 
correlators— in agreement with all previous observations. The new result is that with 
sufficiently hne lattices spacings, i.e., a = 1/8T-1/12T, we hnd a degeneracy of the 
masses from the S/PS correlators and the V/AV correlators. In Figure ^ we demonstrate 
this by showing that the S/PS correlators are coincident with the V/AV correlators up to a 
constant multiplicative factor. 

The local masses plateau for distances z > 1/T and agree with the results of the £t. An 
example is shown in Figure The result of a two mass £t to the PS correlator coincides 
within errors with the plateau in the local masses. Also, the local masses from the parity 
projected components of the V/AV correlators coincide with those from the S/PS 


correlators. Our estimates of the common mass are collected in Table |T1. 

For a = 1/4T our extraction of fips/T is completely consistent with previous estimates; 


for example, a/ipg = 1.144 (4) at T = 2Tc, in perfect agreement with the results of p3 


However, with decreasing lattice spacing we hnd that fips increases rapidly toward the values 
expected for an ideal quark gas, eq. (^). In Figure ^ we show fips/fJ^ppT as a function of 
1/V/ oc a^. Already for a = 1/lOT the screening mass is completely consistent with that in 
an ideal gas. It can be seen that a quadratic extrapolation to the continuum limit at 2Tc is 
consistent with an ideal quark gas. At all three temperatures we see clear degeneracy of all 
the Af~^ screening masses on the larger lattices. The correlators vanish within errors 
even on the coarsest lattice This behaviour is also compatible with expectations from 
an ideal quark gas. 

Interestingly, the measured values of (ijj'ip) differ from the ideal gas value by a factor of 
about two on all the lattices studied. Within statistical errors {'ip'ip) is independent of volume, 
and its ratio to the ideal gas value has small lattice spacing dependence, extrapolating 
smoothly to 1.90 ± 0.04 at 2Tc in the continuum. In view of the chiral Ward identity of 
eq. (P) which relates the pion correlator to the chiral condensate, and the agreement of the 
screening masses with their ideal gas values, we examine the screening correlators for the 
origin of this major departure from ideal gas behaviour. In Figure ^ are shown the PS and 
V screening correlators at T = 2Tc obtained with lattice spacing a = 1/12T. It is clear that 
the correlators are signihcantly different from their ideal gas counterparts, not just in an 
overall factor (which would be sufficient to explain the chiral condensate) but also in shape. 
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FIG. 9: The ratio Xs/Xfft ™ quenched QCD, shown as a function of T. Circles denote results 
for staggered quarks and boxes for Naik quarks. Also shown are the hard thermal loop (HTL) and 
resummed (NL) results from [19|. The bands for the last two are due to uncertainties in Tc/Aj^s- 


V. CONTINUUM PHYSICS 


Our results for the dependence of the continuum extrapolated values of Xs/T"^ oii T/Tc 
are shown in Figure ^ We hnd that the results obtained with staggered and Naik quarks 
are compatible with each other and with the HTL results of at the 2 -ct level. While 
they are not compatible with the resummed results of at the l-cx level, compatibility 
cannot be ruled out at the 3 -ct level for T > 2Tc. There is a suggestion of a drop in 
one approaches T^.. In Table we have listed the continuum extrapolated values of Xs/^^ 
and Xud/T"^ along with the susceptibilities Xs/T"^, Xo/T"^ and Xq/T"^ which are of relevance 
to heavy-ion experiments. 

Previous lattice investigations at cutoff of 1/4T have shown that the effect of unquenching 
is to increase Xs/^^ by less than 


PT| , 0. The computations of indicate that the 
effect of unquenching is to decrease Xa/^^ by less than 5%. Thus the effect of including 
light sea quarks could change the results in Table Q either way within a band of 3-5%. 
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T/Tc 

X3/r' 

XudlT'^ 

XslT^ 

XolT^ 

xqIt'^ 

XsiT) 

p(A++)/27rT 

1.5 

0.84 (2) 

(-2 ±3) X 10“® 

0.53 (1) 

0.43 (1) 

0.99 (2) 

0.63 (3) 

0.96 (1) 

2.0 

0.89 (2) 

(-4 ±4) X 10“® 

0.71 (2) 

0.47 (1) 

1.07 (2) 

0.80 (3) 

1.01 ( 1 ) 

3.0 

0.90 (3) 

(2 ± 2) X 10“® 

0.84 (3) 

0.49 (1) 

1.09 (3) 

0.93 (4) 

1.00 ( 1 ) 


TABLE IV: Results for the continuum limit of quark number susceptibilities, the Wroblewski 
parameter, As(T), expected of a system at thermal equilibrium at temperature T, and the 
screening mass in quenched QCD. We have taken nis/Tc = 1, as appropriate to full QCD. All 
the values reported in this table have been obtained by continuum extrapolation of our staggered 
quark results. Those from Naik quarks are compatible within 1-a limits, except for ^3 at T = l.STc 
which agrees within 95% confidence limits. 


The Wroblewski parameter |^, A^, measures the ratio of primary produced strange to 
light quark pairs. It has been argued [|^] that, for a system undergoing freeze out from 
thermal and chemical equilibrium at temperature T, A^ is just the ratio ‘IXs/iXu + Xd) 
computed at T. We have listed the continuum values of As(T) obtained from our lattice 


simulations in Table |I^. Note especially that As(T) for the staggered and Naik quarks are 
within l-(j of each other. Various extrapolations of Xu = Xd and Xs yield \siTc) ~ 0.4-0.5. 
It seems plausible that sXT k. Tc, Xg dips toward the value obtained from analysis of RHIC 


data, i.e., 0.47 ± 0.04 |34 


Table also contains our determination of the common screening mass in the sector 
for light quarks. At temperatures of 2 Tc and higher this is completely compatible with a 
computation in an ideal gas. Similarly, Xud/T'^ is also compatible with expectations from 
an ideal gas. It has been pointed out before that Xud/T"^ is signihcantly smaller than might 
be expected in the leading perturbative treatment of the interactions. In future it would be 
interesting to investigate perturbative corrections to the screening masses and check whether 
a similar puzzle exists also for them, or whether interaction effects indeed are smaller than 
the error bars on these quantities. Departures from ideal gas behaviour, by a factor of almost 
two, has been observed in the chiral condensate, and hence in the PS susceptibility. It would 
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be useful to have perturbative predictions for these quantities in the future. 
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